Search for: All records

Abstract The resolution and accuracy of single-molecule localization microscopes (SMLMs) are routinely benchmarked using simulated data, calibration rulers, or comparisons to secondary imaging modalities. However, these methods cannot quantify the nanoscale accuracy of an arbitrary SMLM dataset. Here, we show that by computing localization stability under a well-chosen perturbation with accurate knowledge of the imaging system, we can robustly measure the confidence of individual localizations without ground-truth knowledge of the sample. We demonstrate that our method, termed Wasserstein-induced flux (WIF), measures the accuracy of various reconstruction algorithms directly on experimental 2D and 3D data of microtubules and amyloid fibrils. We further show that WIF confidences can be used to evaluate the mismatch between computational models and imaging data, enhance the accuracy and resolution of reconstructed structures, and discover hidden molecular heterogeneities. As a computational methodology, WIF is broadly applicable to any SMLM dataset, imaging system, and localization algorithm. 

We present a computational method, termed Wasserstein-induced flux (WIF), to robustly quantify the accuracy of individual localizations within a single-molecule localization microscopy (SMLM) dataset without ground- truth knowledge of the sample. WIF relies on the observation that accurate localizations are stable with respect to an arbitrary computational perturbation. Inspired by optimal transport theory, we measure the stability of individual localizations and develop an efficient optimization algorithm to compute WIF. We demonstrate the advantage of WIF in accurately quantifying imaging artifacts in high-density reconstruction of a tubulin network. WIF represents an advance in quantifying systematic errors with unknown and complex distributions, which could improve a variety of downstream quantitative analyses that rely upon accurate and precise imaging. Furthermore, thanks to its formulation as layers of simple analytical operations, WIF can be used as a loss function for optimizing various computational imaging models and algorithms even without training data. 

Modulating the polarization of excitation light, resolving the polarization of emitted fluorescence, and point spread function (PSF) engineering have been widely leveraged for measuring the orientation of single molecules. Typically, the performance of these techniques is optimized and quantified using the Cramér-Rao bound (CRB), which describes the best possible measurement variance of an unbiased estimator. However, CRB is a local measure and requires exhaustive sampling across the measurement space to fully characterize measurement precision. We develop a global variance upper bound (VUB) for fast quantification and comparison of orientation measurement techniques. Our VUB tightly bounds the diagonal elements of the CRB matrix from above; VUB overestimates the mean CRB by ~34%. However, compared to directly calculating the mean CRB over orientation space, we are able to calculate VUB ~1000 times faster. 

Simultaneous measurements of single-molecule positions and orientations provide critical insight into a variety of biological and chemical processes. Various engineered point spread functions (PSFs) have been introduced for measuring the orientation and rotational diffusion of dipole-like emitters, but the widely used Cramér-Rao bound (CRB) only evaluates performance for one specific orientation at a time. Here, we report a performance metric, termed variance upper bound (VUB), that yields a global maximum CRB for all possible molecular orientations, thereby enabling the measurement performance of any PSF to be computed efficiently ( $\sim <\#comment/>1000\times <\#comment/>$faster than calculating average CRB). Our VUB reveals that the simple polarized standard PSF provides robust and precise orientation measurements if emitters are near a refractive index interface. Using this PSF, we measure the orientations and positions of Nile red (NR) molecules transiently bound to amyloid aggregates. Our super-resolved images reveal the main binding mode of NR on amyloid fiber surfaces, as well as structural heterogeneities along amyloid fibrillar networks, that cannot be resolved by single-molecule localization alone.